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Abstract : 

Specific heat is calculated using Tsallis Statistics. It is observed that it 
is possible to explain some low temperature specific heat properties of glasses 
using non-extensive approach. A similarity between temperature dependence of 
non-extensive specific heat and fractal specific heat is also discussed. 
Motivation: 

Non-extensive statistics is being increasingly used to explain anomalous be- 
haviour observed in the properties of various physical system. Tsallis statistics 
has been used to study physical systems / phenomena which include turbulence 
in plasma [1] , Cosmic ray background radiation [2] , self gravitating systems [3] , 
econo-physics[4], electron -positron annihilation [5], classical and quantum chaos 
[6], linear response theory [7], Levy type anomalous super diffusion [8], ther- 
malization of electron - phonon systems [9] , low dimensional dissipative systems 
[10] etc. It has been shown that non-extensive features get manifested in those 
systems which have long range forces, long memory effects or in those systems 
which evolve in (non Euclidean like space-time) fractal space time [11 and ref- 
erence therein]. Apart from other applications it has been suggested [11] that 
non-extensive statistics can be applied to complex systems like glassy materi- 
als and fractal / multi-fractal or unconventional structures also. Anomalous 
low temperature specific heat results in glasses have motivated us to use non- 



extensive statistics. We will also compare specific heat using Tsallis statistics 
with specific heat of a fractal. 
Low temperature specific heat: 

Specific heat depends on density of states g(co) 

where x =-j^pf- Above equation has been obtained using Boltzmann Gibbs 
statistical mechanics. Glasses at low temperature (below 2K) show quasi-linear 
behaviour i.e. approximate linear dependence on temperature T. Again glassy 

systems do not follow T 3 dependence even above 2K as in some cases maximum 

c 

in j% is observed. Above properties are universal features of glassy systems 
[12,13]. Another interesting feature is ^f- dependence on T 2 for crystals and 
glassy phase of the same systems [14]. 
Theory : 

Non-extensive statistics is based on two postulates [11 and reference therein]. 
First is the definition of non-extensive entropy 

1-P? 

S q = — f ( 2 ) 
where q is the non-extensive entropic index and Pj are the probabilities of 

the microscopic states with ^2 Pi = 1. 

The second postulate is the definition of energy U q — PfEi, where Ei is 

the energy spectrum. As q — >1, S q — - ^2pilnPi, which is the Boltzmann Gibbs 

Shannon entropy. 

To derive non-extensive form of specific heat we have to use non-extensive 
quantum distribution function of bosons. It is very difficult to derive exact 
analytical expression for non-extensive distribution function. However, there are 
many studies which provide the approximate form of non-extensive distribution 
functions [2,15]. In the present paper we use dilute gas approximation (DGA) 
for boson distribution function. For DGA case, the average occupation number 
is given as [15] 

1 



(l + (q-l)f3(Ei-^ -1) 



(3) 
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Figure 1: Non-extensive specific heat for various values of q; 
q=1.5(a), 1.6(b), 1.7(c), 1.8(d), 1.9(e),2.0(f),2.5(g),3.0(h),3.5(k) 



When dealing with system in contact with the heat bath at the temperature 
/3, we have 

< n q >= r — i 4 

(l + (g- 1)^)3=1-1 

Averaged over non-extensive distribution, expectation value of energy is 

< E q >= hui <n q > (5) 
and total vibrational energy is 



J fiuj < n q > g(uj) (6) 



from which specific heat can be obtained 

C - ™1 (7) 

For the case of Debye approximation g(u) cx (uj) 2 . Equation (7) represents 
specific heat using Tsallis statistics where n q is given by equation (4). 
Results and Discussion : 

We have numerically solved equation (7) using equations (4) and (6), for 
various values of q. It is very clear from figure 1 that quasi-linear behaviour 
or maxima in ^ is not explained by including non-extensive statistics. Vari- 
ous curves in figure 1 have been obtained for fixed value of Debye temperature 
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Figure 2: Fig2a: Sp. heat data corresponding to {B 2 0z)im-x{Na 2 0) x for 
different % of Na 2 [16,17,18] Fig2b: Sp. heat data corresponding to glassy 
alcohols [14,19]. In y-axis ^ is plotted against T 2 in x-axis 

(#d=300 K). The shape and nature of curves in figure 1 is decided by value of 
q. A linear display of specific heat versus T 2 dependence is followed in glasses 
as depicted in figure 2a and figure 2b. In figure 2a we are plotting experimen- 
tal data corresponding to (B 2 3 ) W o- x (Na 2 0) x taken from reference [16,17,18]. 
The data with (+) sign corresponds to % of Na 2 0, (*) corresponds to 1 % of 
Na 2 0, hollow square corresponds to 6 %, filled square corresponds to 16 % and 
crosses (x) corresponds to 25 % of Na 2 respectively. The fitted curves corre- 
spond to different values of q from (+)onwards are q=5.5, 5.0, 4.5, 3.5,3.0,2.9, 
2.7,2.6 respectively. The slope of experimental data in figure 2a is strongly de- 
pendent on presence, absence or excess of Na 2 in glassy system. Non-extensive 
specific heat seems to explain data in figure 2a for different values of q. Again 
in figure 2b specific heat data of various glassy systems show strong dependence 
on slope and can be explained by non-extensive approach. The experimental 
data in fig2b corresponds to glassy alcohols like 2-proponal (+), filled square 
to D-ethonal, cross (x) to D-ethonal (OG, orientational glass), hollow square 
to 1-proponal and star (*) to glycerol respectively. The corresponding fitted 
q-values are q=2. 8, 2. 65, 2. 62, 2. 60 and 2.25 respectively. The experimental data 
has been taken from reference [19,22]. It is interesting to mention here that a 
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Figure 3: Comparison of non-extensive specific heat with fractal specific heat. 
In y-axis ^ is plotted against T 2 in x-axis 

strong slope dependence of specific heat on T 2 has been used to [14] depict two 
different phases of crystal and glasses with two different Dcbye temperatures. 

Low temperature specific heat studies in silica aerogels and other disordered 
systems do not follow Debye model and concepts of fractals have been used to 
explain such results [20] . For a fractal approach density of states g(u>) oc uj d 
where d is spectral dimension. 

We have plotted proportional temperature dependence of specific heat of a 
fractal in figure 3. Curves b,d and e correspond to d=3.8,2.8 and 1.8 whereas 
curves a,b and f correspond to q=1.7,1.65 and 1.6 respectively. It is evident from 
figure 3 that there is lot of similarity between fractal and non-extensive specific 
heat. The similarity between these curves is expected because non-extensive 
statistics corresponds to non-Euclidean space time which is fractal in nature. 
The similarity between fractal and nonextensive results were also reported for 
Lamb Mossbauer factor [21]. It has been shown that magnetic properties of 
manganites [23] and their corresponding phase diagrams [24] can be explained 
in the framework of Tsallis statistics. 

Some of the anomalous properties have been explained by using soft poten- 
tial model which is an extension of tunnelling model in which group of atoms 
can tunnel between same energy configurations [14]. Recently it has been re- 



ported that quasi-linear behaviour or maximum in specific heat data of various 
incommensurate phases can be explained [25] by considering effect of phason 
damping. 
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